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Jean-Pierre Eckmann; Guido Schneider 

Nonlinear Stability of Bifurcating Front Solutions for the Taylor- 
Couette Problem 

We consider the Taylor- Couette problem in an infinitely extended cylindrical domain. There exist modulated 
front solutions which describe the spreading of the stable Taylor vortices into the region of the unstable Cou- 
', ette flow. These transient solutions have the form of a front-like envelope advancing in the laboratory frame 

' and leaving behind the stationary, spatially periodic Taylor vortices. We prove the nonlinear stability of these 

^ solutions with respect to small spatially localized perturbations. 
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1 Introduction 



Ludwig Prandtl's boundary layer theory of 1904 was a breakthrough in modern fluid mechanics. It allows to 
distinguish two domains for small viscosity flows. In a narrow neighborhood of the boundaries viscous effects play 
Ph a role while outside this boundary layer the flow can be considered as a potential flow. The Taylor-Couette problem 
—< ' we consider here is a hydrodynamic problem where the onset of instability is due to the motion of the boundaries. 
^ i ' While Prandtl's ideas are still relevant, they do not apply as such in our case because the infinite length of the 

boundaries will produce an infinite amount of energy. 
^ I The Taylor-Couette problem consists in describing the fiow of a viscous, incompressible fluid in the domain 

^ ' between two (counter-) rotating cylinders. More precisely, we are interested in the velocity field of a viscous 
incompressible fiuid filling the domain = R x S between two concentric rotating infinite cylinders, where E C 

■ denotes the bounded cross section. The flow between the rotating cylinders is described by the Navier-Stokes 
^ equations in Q with no-slip boundary conditions. For simplicity we assume throughout this paper the outer 

■ cylinder to be at rest. 

' When the rotational velocity of the inner cylinder is small, a stationary flow appears, called the Couette 

flow C/cou with associated pressure flcld pcou- It is homogeneous along the cylinders and the streamlines are given 
by concentric circles. For small Reynolds number TZ (which is proportional to the rotational velocity of the inner 
cylinder) these solutions are exponentially stable. 

A perturbation {U,p) of the Couette flow satisfies the Navier-Stokes equations 



dtU = A^7-7^[(^7cou•v)^7 + (c/•v)^7cou + (c/•v){/]-Vp, (i) 

V • t/ = , (2) 
with Dirichlet boundary conditions 

C/kxas = . (3) 
This problem has a unique solution U and Vp if we add the vanishing mean flux condition. 



[U(^)h{x) = -p^ / C/(x)(a;, z)dz = , for aU X e M , (4) 

1^1 -'zGS 



where C/(x) stands for the velocity component along the infinite x-direction. See g. 

The trivial branch of solutions of the system (0)-(0) is U{x, z) = 0, which corresponds to the Couette flow. 
It becomes unstable when the Reynolds number TZ goes beyond a certain threshold TZc. We define e > by 
e'^ =TZ — TZc and assume throughout £ <C 1. 

For e > a family of small spatially periodic equilibria (Taylor vortices) bifurcates from the trivial solution 
U = 0: 



C/TV,[g,e,^](a;, z) ^sA^l- c,q^ [/^^ (z)e*('=^+-'?)(-+'') + c.c. + 0{e^), qe[0,l/^), deR. (5) 

In (||), ci is a positive constant, /cc > is a critical wavenumber, and Uk^{z) £ C"^. We will henceforth assume 
•d = and omit it from the parameters of Utv- Such solutions have been constructed in pq, p53, Eq]. It turned out 
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that the hnearization around the stationary solutions f/Tv,[(j,e] possesses continuous spectrum at least up to the 
imaginary axis. These solutions are linearly stable for ciq^ < i and so-called sideband- or Eckhaus unstable for 
ciq2 > i. See or §. 

In actual experiments, the transition from the Couette flow to the Taylor vortices usually begins from per- 
turbations at the ends of the (finite) cylinders. A first vortex forms near the end of the cylinder, and from this 
seed new vortices form away from the ends, and an advancing front seems to move through the cylinder, leaving 
stationary, spatially periodic vortices behind. Finally the cylinder is filled with Taylor vortices. 

These advancing fronts cannot just be described as U(x,t) = Uf{x — ct), because of the pattern they leave in 
the laboratory frame. Rather they are of the form described in 



Theorem 1. J^/ . For each q with ciq^ < ^ the following holds. For sufficiently small e > and for c > 
the system ©-(I) possesses modulated front solutions 



with 



U{x,Z,t) = Un-,f,[q,e,c]ix - ct,X,z) , 



Uml,[q,e,c] i^i X, z) — J7nif,[g,e,c] + 



27r 



eq 



and the boundary conditions 



lim U^i,[q,e^c]{^,x,z) = C/Tv,[g,e](a;,z) 



lim ?7mf,[g,£,c](C,2:,z) 
4 — *-oo 



= 



the Taylor vortices , 
the Couette flow . 




eriveippe(i') / X\ 

/ \ V eirvelope(t) 



\ / \>CAnf,[g,£,c](a;- ci,a;,z) 

pml,[q,sk(i^ ~ Ct',X,z) 



Figure 1: Sketch of a modulated front, advancing from left to right at times t' < t, and the advancing envelopes. 

In the present paper we prove the nonlinear stability — under spatially localized perturbations — of these bifurcated 
front solutions f/inf,[g.e,c] • In particular, this means that a small perturbation cannot destroy the growing regular 
pattern in the bulk, nor make it move. This is not an obvious result since there could be a soft mode in the 
system, due to translation invariance. But, although the Taylor vortices are only diffusively stable, their positions 
are pinned. The precise statement of the result can be found in Theorem ^. 

Our proof is an adaptation of our previous paper JTol where we dealt with the stability of front solutions for 
the Swift-Hohenberg equation 



dtu = -{I + d'J^u + 



- e^u - 



u{x,t) e 



X e 



t > 



connecting the stable periodic equilibria 



= £ 



L_i^e^(l+e9)-+C.C. + 0(£2) 



for Aq^ < 1/3 with the unstable trivial solution u = 0. 

The stability proof is based on the following idea. Consider a spatially localized perturbation ahead of the 
front. In this unstable region any perturbation will grow exponentially because is an unstable solution. But, after 
a finite time the front will have reached the perturbation, or equivalently, the perturbation has entered the stable 
region behind the front. This region is diffusively stable and thus perturbations decay as 1/Vt. Therefore, one 
expects the front to be stable under sufficiently small spatially localized perturbations which decay exponentially 
with the distance ahead of the front. In particular, the front is stable for sufficiently small perturbations with 
compact support. 
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In order to make these ideas clearer, we describe first two easier model problems. In Section ^ we explain for 
a very simplified model how weighted norms can be used to stabilize the a priori unstable region ahead of the front. 
In this very simplified model the equilibrium in the bulk will be assumed to be stable with some exponential rate. In 
Section ^ we give up the exponential stability of the equilibrium and replace it by a diffusively stable equilibrium, 
i.e., the continuous spectrum of the linearization now reaches the imaginary axis and is no longer separated from 
it. As an example we consider the amplitude equation associated with the Taylor-Couette problem, namely the 
Ginzburg-Landau equation. For this equation it will be more or less obvious that spatially localized perturbations 
of the stable equilibrium decay as a solution of a linear diffusion equation. For the Taylor-Couette problem, this is 
no longer obvious for the solutions we are interested in and so the main difficulty for the fronts C/mf is to recover 
the diffusive behavior behind the front. 

Remark 1. Since the internal wavelength of the spatially periodic flow is in general much smaller than the 
length of the cylinders and due to the fact that the ends of the cylinders in a laboratory experiment do not influence 
the form of the pattern in the interior, we consider the Taylor-Couette problem with cylinders of infinite length. 
As explained, up to a certain extent this idealization is a good description of Nature. 

Remark 2. The Taylor-Couette problem is one of several examples where our theory applies. Other exam- 
ples are Benard's problem and reaction-diffusion systems. Loosely speaking our theory applies whenever a spatially 
homogeneous state becomes unstable via a supercritical stationary bifurcation at a non-zero Fourier wavenumber 
in a translationally invariant system. 

Remark 3. The existence proof for modulated fronts U{x,z,t) = C^nf, [(?,£, c] (2^ ~ ct,x,z) connecting the 
stable Taylor vortices with the Couette flow is based on center manifold theory for elliptic problems on unbounded 
cylindrical domains p^ . It is complicated here by an infinite number of eigenvalues on the imaginary axis for 
TZ = TZc- The related result for the Swift-Hohenberg equation can be found in pjf . 



2 Basic ideas 

The ordinary differential equation dtu — u with u{t) e M is used to model the exponential growth of a certain 
quantity u, for instance a chemical concentration or a population of bacteria. Usually nonlinear saturation effects 
give some bound on the possible size of u. This can be taken into account by adding a nonlinear term, for example 
— M^, to the model, leading to 

dtU = u — . (6) 

This nonlinear model possesses in addition to the unstable equilibrium u = also a stable equilibrium u = 1, and 
so every solution with initial condition in the interval (0, 1) converges towards u = 1 for t — s- 00. This model is 
refined further by taking into account the spatial structure of the solutions: Chemical reactions or the growth of a 
population in general do not happen in a spatially homogeneous manner and so the above model is extended by 
adding a diffusion term d^u to the system, i.e., 

dtu — d^u + u — V? , (7) 

with u{x, t) G R. For solutions with small positive spatially localized initial conditions the dynamics is very simple. 
Solutions with such initial conditions develop in a universal way. By the reaction term u — v? at fixed a: G R the 
solutions are drawn to the stable equilibrium u = 1. Then, by the diffusion term the part of the solution which 
is close to the stable phase u = 1 starts to spread into the region of the unstable phase u = forming two fronts 
with velocities c = ±2, one to left and one to the right. 

Such front solutions u{x,t) = Wf [c](a; — ct) = Wf,[c](C) connecting the stable phase u = I for ^ — > —00 with the 
unstable phase u = for ^ ^ 00 can be found by analyzing the second order ordinary differential equation satisfied 
hy V = ffjc] in the {v,v')-pha,se plane. The partial differential equation (|^) can be considered in the moving frame 
£^ — X — ct, where it takes the form 

dtv = d^v + cd^v + V — v'^ . (8) 

A front Wf [j,] then satisfies v'/^c] +ci'£ [^j +^f,[c] [c] ~ ^- '^^^ found by showing that there exists a heteroclinic 
connection between (w, v') — (1, 0) and {v, v') = (0, 0), where the latter fixed point has real eigenvalues when |c| > 2. 
As a consequence, when |c| > 2 the fronts satisfy Wf,[c](0 G (0; !)■ 
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Thinking of a chemical reaction front it is obvious that the stabiHty of these fronts is a relevant question. We 
recall that an equilibrium y = of an ordinary differential equation y — Ay + 0{y^) is called (Lyapunov) stable in 
if for all e > there exists a (5 > such that ||y(0)||Rd < 5 implies ||?/(t)||Rti < £ for all t > 0. 

A good starting point is to consider stability under perturbations in the space of uniformly continuous 
bounded functions equipped with the norm ||u||co = sup^g^ Since the fixed point m = is already unstable 

for the ordinary differential equation (||) it is also unstable under the dynamics of the partial differential equation 
(^. On the other hand, perturbations of m = 1 in can be bounded by the maximum principle using spatially 
constant functions, and so the stability of this equilibrium follows again from the stability of u = 1 in the ordinary 
differential equation (^). 

A drawback of this topology is that the fronts I'fjc] s-re actually unstable with respect to perturbations in C^. 
In order to see this consider the initial condition 0) = max(t)(^), 5). Then, as we see in Fig.^, 

sup 0) - < <5, but lim sup t) - = 1 . 



However, we shall now argue that these fronts are stable with respect to spatially localized perturbations. Suppose 
that a perturbation is localized near ^ = > 0. Then the perturbation has only a finite time t = ^o/c to grow 
exponentially before it is hit by the bulk of the front, where for (0) perturbations are damped exponentially. 




exponential decay 




X 

exponential growth 

Figure 3: Stability mechanism for fronts of (^. 

Next we explain how this heuristic argument can be made rigorous with the help of exponentially weighted 
norms. In order to do so we consider the stability of the zero solution u = for the equation dtu = 9|u + cd^^u + u. 
Let u'(^) = u{£^)e^^ with a /3 > 0. Then the function w satisfies 

dtw = d^w + (c - 2(3)d^w + (/3^ - + l)w . 

If — /5c + 1 < (which we assume henceforth) the maximum principle implies the stability of w = in C^. 

However, this reasoning does not carry over to the nonlinear problem since the product of two functions 
Uj{x) — Wj{x)e~^^ , {j = 1, 2) with Wj G leads to uiU2e^^ = wiW2e~^^^e^^ , which is in general not bounded for 
X —^ —oo. Therefore, the exponential stability of w = cannot be used directly. 

Only an interplay of both norms, i.e., the consideration of u and w simultaneously, and the presence of a 
front will suffice to show the stability. We consider thus a front Wf traveling with speed c, which is a solution of 
(resp. (||)). The deviations r = u — Uf [cj and w{^) = r(^)e^^ from the front Uf^jc] satisfy 

dtr = d^r + cd^r + r — 2vf, [c]i^ — , 

dtW = a|w + (c - 2P)d^w + {(P ~ I3c+ l)w - 2vf ^[c]W - rw . 
At first sight this seems to bring no advantage, but the term — 2uf [^jr in the equation for r can be expressed as 

-2vi^[c]r = -2r + 2(1 - ^;t,[c])f = -2r + 2(1 - V[^c])e''^^e'^^r = -2r + 2(1 - vi^[c])e~^^w . 
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Thus, the equation for r is transformed to: 

dtr = 9|r + cd^r — r + 2(1 — vi^[c])e~^^w — . 

The term 2(1 — 'i'f.[c])e^'^^ is hounded when — /3c + 1 < 0, since 1 — Wf,[c](C) = ©(e'^'^l) as ^ — > -co with 7 > 
solving 7^ — 7c — 1 = 0. Because of the term — r in the equation for r and because we have — 2uf [c] < in the 
equation for w the hncarization of the system (r, w) around (r, w) — (0, 0) is now exponentially stable and so this 
obviously also holds for the nonlinear system. 

Theorem 2. Assume /3 and c satisfy (3^ — f3c + 1 < 0. Then for all e > there exist C\S > such that 
sup(|r(e,i)| + \w{^,t)\) < Ce('°"''(-i''''-'3^+i)+")* , 

provided 

sup|r(e,0)(l + e'3«)| < S . 

The discussion of this equation goes back to ^ . The above result can be found for instance in ||3^ . 

Remark 4. The discussion of the marginally stable case — c/3 + 1 = is much more delicate, but covers 
the important case of minimal velocity |c| ~ 2. For these solutions global stability results have been obtained in 
. The exact asymptotic decay of the perturbations in the sense of the next section has been investigated in j2^ . 



3 The Ginzburg-Landau equation 

In this section we consider a model which comes even closer to the Taylor- Couette problem, namely the Ginzburg- 
Landau equation: 

dtA = dlA + A-\AfA, (9) 
with A = A{x, i) e C, t > 0, and a; e R. This equation has a two-parameter family of nontrivial equilibria 

^cq,M](a;) = v/r^e'^'"+''^ • 

In his pioneering work, Eckhaus [Q showed that the equilibria with < 1/3 are linearly stable while those with 
> 1/3 are linearly unstable. In contrast to the equilibrium w = 1 of the previous problem the stable equilibria 
^cq,[q,^] s-re no longer exponentially stable. Indeed, the linearization always possesses continuous spectrum (—00, 0] 
up to the imaginary axis, i.e., at best polynomial decay rates can be expected. 

The following discussion can be carried out for all < 1/3. Since our aim is to explain the basic ideas, we 
will restrict ourselves to the stability of the front A{x, t) = B{x — ct) connecting the stable equilibrium A(,q jq gj = 1 
with the unstable equilibrium A = Q which corresponds to the case q — Q, ^ — Q. 

The real-valued front B = i?f.[c](C) satisfies the ordinary differential equation B" + cB' + B — B^ — 0. It 
can be found by analyzing this second order ordinary differential equation in the {B, _B')-phase plane. There 
exists a heteroclinic connection between {B, B') = (1, 0) and {B, B') — (0, 0), where the latter fixed point has real 
eigenvalues for |c| > 2. As a consequence, when |c| > 2 the fronts satisfy -B(^) € (0, 1). 

We introduce the deviation v from the front by writing A{x, t) = i3f,[c] [x — ct) + v{x, t), and if A satisfies (^) 
then V solves the equation 

dtv = dlv + V ^ Bl^^-^{2v + v) - v^B^,] - 2\vfBf^[,] - v\v\^ . (10) 

Again, the stability of the front will follow only by an interplay of the usual norm and the spatially exponentially 
weighted norm. Therefore, we introduce w{^) = + ct)e^^ and the equation for w is 

dtw = dlw + (c - 2(3)d^w + {P^ -cf3+ l)w 

-Bf [^](2w + w) - B{^[c]vw - 2B{^[c]Vw - \vfw , 

where ^ — x — ct. When — c/3 + 1 ~ —2j < 0, the maximum principle applied to the system written in polar 
coordinates w = \w\e^'^ implies the exponential stability of w = and we get a bound w = 0{exp(^jt)) as t — s- cx) 
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provided v stays bounded. To exploit the effect of the bulk, we write the difference between i?f jc] and the stable 
equilibrium A(,q.[o.o] = 1 as g = Bf [j,] — 1. Using g, ( |lo|) can be written as 

dtv = dlv + v-{2v + v)-v'^ ~2\v\^ -v\v\^ 
^{2g + g''){2v + v)~v^g~2\v\'g . 

All terms in the second line vanish at some exponential rate as x ~ ct —oo. We use this as above by expressing 
V in terms of w, for instance, gv — {e~^^g){e^^v) — {e~^^g)w, and so 

dtv = dlv + v- {2v + v)~v'^ - 2\vf - v\vf + 0{w) . 

Writing A in polar coordinates A — {1 + r)e*''' we then find 

dtr = dlr ~2r- (d^ipfil + r) -Sr"^ -r^ +0{w) , (11) 

dtV = dxV + ^ — r- +0{w) . (12) 

1 + r 

Linearizing this system, ignoring the terms 0{w) which decay at some exponential rate, gives 

dtr = d^r — 2r and dtf — d^^p . 

This shows that the variable r decays with some exponential rate. However, the variable ip shows some diffusive 
behavior. Therefore, the variable r can be expressed asymptotically by (p, that is, r = —(dx^)'^ /2 + h.o.t. . Inserting 
this into the equation for ip yields 

dt(p = dlip - {dx(p)dx{dx(ff + h.o.t. . 

We explain now the arguments which are needed to analyze this equation. In order to do so we go back to the 
linearized system. In Fourier space ip satisfies dt<fi{k,t) = —k'^(f{k,t), that is (p{k,t) — exp(— fc^i)(^(fc, 0). Thus, the 
Fourier modes of (p are concentrated around fc = as time evolves and so 

(p{k/Vt,t) = eyip{-k^)>p{k/Vi,0) = exp(-fc2)^(0, 0) + C'(l/\/t) , 

for sufficiently smooth (p{k,0). Since smoothness in Fourier space corresponds to decay rates for — *■ oo this 
explains the well known fact that solutions to spatially localized initial conditions decay to zero in a universal 
manner, 

^(a:,i) = ^exp(^) + 0(t-i), 
with some constant (/s, e R depending only on the initial conditions. As a consequence, we have 

but under these asymptotics. 

Therefore, all nonlinear terms in the equations for (r, tp) vanish much faster than those which are part of the linear 
diffusion equation. In general, nonlinear terms (^^^ (9a;y)P^ (i9^(/?)^^ are therefore called irrelevant if pi+2p2+3p3 > 3. 

With the help of renormalization theory these heuristic arguments have been made rigorous and it has been 
shown that the nonlinear system possesses the same asymptotics |l^, |^ . Thus spatially localized perturbations of 
the equilibrium ^oq,[o,o] = 1 vanish asymptotically like solutions of a linear diffusion problem. 

Notation. Let ||u||l2 = / \u{x)\'^dx. We define the usual Sobolev norm by ||m||h2 = II^Hl^ + ||92:u||l2 + 
||9^u||l2. Moreover, we set ||u||h2 = ||up||h2, where p{x) = (1 + x^y/^. 

Theorem 3. There exist positive constants 6, and C , such that for all vq with \\vo\\-^2 < S the following 
holds. Let A = A(,q_[o,o] + v be the solution of the Ginzburg- Landau equation (|^) with initial condition x i— > 
^cq,lo,o]{x) + vq{x) = l + t;o(a;)- Then v exists for all times t > and there exists a constant € M depending only 
on the initial condition vq such that for all t > 1 one has 

sup\v{x,t) - ^exp(— ^)| < ^ . 
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Furthermore, spatially localized perturbations of the front A{x, t) = Bf\^^-^{x—ct) vanish asymptotically as a solution 
of a linear diffusion equation |l^ . 

Theorem 4. Let Sf,[c] be a front with velocity c > 2. Let f3 > be such that — c/? + 1 = —27 < 0. 
Then there exist positive constants S and C such that for all vo with ||x > vo{x){l + e'^^)||H| < <^ the following 
holds. LetA{x,t) = B{ i^c^{x — ct) + v(x,t) be the solution of the Ginzburg-Landau equation with initial condition 
Aq(x) = i?f + Vq{x). Then v exists for all t > and there is a constant depending only on vq such that for 
all t > 1 one has 

sup|w(x,t) - ^exp(-^)| < ^ , 

x£R \Jt 4i t 

and 

sup|f(C + ct,i)e'^«| < Ce""^* . 

The case of minimal velocity c = 2, which is more delicate, has been handled in ^ with norms which take more 
details of -Bf [^j into account. 

Remark 5. These fronts are closely related to the modulated fronts of the Taylor-Couette problem. Indeed, 
using the multiple scaling ansatz ^ , 

U ^'\l)e,A + 0(e^), where = £ ^(ea;,e^t) C/fc^ (x)e*'''^^ + c.c. , 

one sees that the Taylor-Couette problem is approximated by the Ginzburg-Landau equation for the A. Mathemat- 
ical theorems which show exact relations between the Ginzburg-Landau equation and the Taylor-Couette system 
can be found in §, ||, |2|, ||, |3|, 0. 



4 The Taylor-Couette problem 

Li line with the discussion of the previous sections the stabihty proof for the modulated fronts Uyai\q,z,c\ (a; — ct, x, z) 
mainly consists of two parts: i) The introduction of an exponential weight to stabilize the unstable part ahead of 
the front and ii) A stability analysis of the stationary solution V'yv \q.e\ behind the front. As above we will restrict 
ourselves in this paper to the case g = which lies in the Eckhaus-stable region. 

Note that for this problem the stability analysis of the equilibrium is by no means obvious, but can be found 
in ]4l| ]. The equilibrium C/TV,[o,e] of possesses continuous spectrum up to the imaginary axis just like 

^cq,[o,o] of Eq.(|^). Therefore, one can expect at best polynomial decay in time. Furthermore, it has been shown that 
the linearly stable Taylor vortices are also nonlinearly stable with respect to small spatially localized perturbations 
and that the perturbations vanish in a universal way like a solution of a linear diffusion equation. 

Theorem 5. Let £ > sufficiently small. Then there exist S > 0, and constants Ci > 0, C2 such that 

for all Vo with ||Vo||fj| < <5 the following holds. Let U — ?/TV,[o,e] + V he the solution of the system (|l|)~(^) with 
initial condition {x,z) t-^ t^TV.[o.e] (^^j ^) + Vaix,z). Then V exists for all times t>0, and there exists a constant 
S R depending only on the initial condition Vq such that for all t > 1 one has 

1 —x^ C2 

sup \V{x,z,t) -=V;exp(— — )a^C/TV,[o,e](a;,z)| < . 

{x,z)eKxT. yt 40it 

The corresponding result for the Swift-Hohenberg equation can be found in |l^ . 

After the stability analysis of the part behind the front it remains to stabilize the zero solution ahead of the 
front. We again use the interplay of two norms. 

We introduce W by z, t) = -I- ci, z, t)e^^. If we insert this into the linearization around the Couette 
flow we obtain an estimate for the variable W , namely 

\\W{;;t)\\n^ < Ce~'^^''0^'^'\\W{;;O)\\H^ ■ 

For given e > 0, sufficiently small, there is a minimal velocity cq > such that for all c > cq there are f3i and (3u 
such that for all /3 £ (/?;,/?„) one has p(c, /3, e) > 0. Given such an e > and a c > co we choose now a f3 such that 

p{c,P,e) = 27 > . 



Then we have: 
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Theorem 6. Let e > sufficiently small. Then there exist a 6 > 0, and constants Ci > 0, C2 such that for 
all Vq with II (a;, z) i~> Vo{x, z)(l + e^'^)\\ii2 < S the following holds. Let U{x, z,t) = U^{^[o,e,c]ix — ct, x, z) + V{x, z,t) 
be a solution of the system with initial condition [x, z) i— > C/mf^[o,e,c] (a^i x, z) + Vo{x, z). Then V exists for 

all times t > 0, and there exists a constant 14 G M depending only on the initial condition Vq such that for allt>l 
one has 

1 —x'^ C2 

sup \V{x,z,t)-—^V^eiiY>{jTr-;)dxUTY^iQ^e]{x,z)\<-—, (13) 

and 

sup \V{^ + ct,z,t)e^^ <C2e^^* . 

(5,z)eRxS 

Sketch of proof. We proceed as in the previous sections. First we write as a dynamical system in the 

space of divergence free vector fields |^. Then we consider the system for V and the exponentially weighted 
variable W. The proof is exactly the same as in |jl^ given for the fronts in the Swift-Hohenberg equation. The 
analysis behind the front and the complete functional analytic set-up can be found in |^. Finally, in |4^, one 
finds the estimate showing that the Taylor-Couette problem is approximated by the Ginzburg-Landau equation. 
This is used ahead of the front. □ 



Remark 6. One can interpret (y_3|) as saying that the perturbations decay faster near the extrema of Utv- 
This can also be understood in the sense that at large amplitude of the underlying periodic pattern (i.e., where 
dxUTv{x, z) = 0), the restoring force is stronger. 

Remark 7. The choice of a sufficiently small e > actually allows to prove the stability of all fronts which 
are predicted to be stable by the associated amplitude equation (the Ginzburg-Landau equation) since 

lim e"^e(ec_B,e/3A,e) = e^(c_B,/?A) , 

where ^^(cb, /3a) is the stability condition of the associated amplitude equation. 

Remark 8. Note that the Taylor-Couette problem is structured by the periodic background of the stationary 
solution. In such a setting it is natural to study perturbations (which need not be periodic) in the so-called Bloch 
wave representation (see e.g. Q). This means that one writes 

.1/2 

u{x) ^ / e''^'=u{£,x) , 
J-1/2 



where 



with u the Fourier transform of u. Since the linearized problem for the Taylor vortices is periodic, the corresponding 
linear operator is diagonal in the £ and thus the Bloch representation simplifies the analysis. Nonlinear terms can 
be written in Bloch space with the help of convolution integrals. Decay for \x\ — s- 00 corresponds in Fourier and 
Bloch space to smoothness with respect to the Fourier wave number k and Bloch wave number £, respectively. As 
a consequence, the important relation dx- = 0{t~^^'^) for the irrelevance of terms corresponds in Bloch space to 

f =0(i-l/2). 

Remark 9. The diffusive behavior of solutions close to the Taylor vortices can be understood as follows: The 
eigenfunctions of the linearization of the Taylor vortices are given by Bloch waves e'^^^Vn{£, x, z) with u„(£, x, z) = 
Vn{(,x + 27r/fcc, z) having the same periodicity as the Taylor vortices. For each fixed Bloch wave number £ there 
is discrete spectrum, i.e., n € N. Therefore, there are smooth curves fj.ni£) of eigenvalues over the Bloch wave 
numbers. 

There is one curve, ni, coming up to zero. This curve has approximately the form of a parabola, i.e., 
fii{£) = —cil"^ + 0{l^) with ci > in the Eckhaus stable region. Thus, the solution -Oi rescaled in Bloch space 
satisfies 

lim vi{e/Vi,t) = u,e~'=i^' , 

t — *oo 

for smooth initial conditions in Bloch space, and shows the same asymptotic behavior as the linear diffusion 
equation dtVi = Ci&^vi for spatially localized initial conditions. 
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Figure 4: Spectrum of the linearization around the Taylor vortices as a function over the Bloch wave numbers. 

Remark 10. As for the case of the Ginzburg-Landau equation we have, for the perturbations of the Taylor 
vortices, asymptotically irrelevant nonlinearities. This is not obvious and has been proved in |4l| ] by deriving 
an effective equation for the variable corresponding to the curve of critical eigenvalues fii and proving that the 
coefficients in front of the quadratic and cubic nonlinear terms vanish up to a certain order as the Bloch wave 
number £ goes to zero. These facts are a reflection of the translation invariance and of the existence of a circle of 
fixed points for the stationary problem. 

Remark 11. The nonlinear stability of so-called critical fronts (moving at the minimal possible speed for 
which they are linearly stable) remains open. See for the linear stability analysis of the fronts in the Swift- 
Hohenberg equation. Achieving this aim seems to be a necessary step in solving the long-standing problem of 
"front selection" [T^ , in a case where the maximum principle |l] is not available. 
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